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Abstract 



It is proved that the roots of combinations of matrix polynomials with real 
^^ ■ roots can be recast as eigenvalues of combinations of real symmetric matrices, 

^^ . under certain hypotheses. The proof is based on recent solution of the Lax con- 

jecture. Several applications and corollaries, in particular concerning hyperbolic 
matrix polynomials, are presented. 
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^ • 1 Main result 

—I. I A polynomial is called hyperbolic if all its roots are real. It is a classical well studied 

O I class of polynomials (see, e.g., JE])- There are at least two useful ways to extend 

^ ■ this notion to polynomials with complex n x n matrix coefficients, in short, matrix 

d • polynomials. Thus, a monic (i.e., with leading coefficient /„, the nxn identity matrix) 

matrix polynomial L{z) of degree i is said to be hyperbolic, if for every nonzero x G C^, 

the n-dimensional vector space of columns with complex components, the polynomial 

equation 

'C- {L{z)x,x)=0 (1.1) 

has i real roots (counted with multiplicities). We denote here by (■, ■) the standard 
inner product in C". An nxn monic matrix polynomial L{z) of degree i will be called 
weakly hyperbolic if detL{z) = has ni real roots (multiplicities counted). Note that our 
terminology differs slightly from the terminology in some sources (for example, [Ti]). 
Clearly, every hyperbolic matrix polynomial is weakly hyperbolic, and the coefficients 
of every hyperbolic matrix polynomial are Hermitian matrices. See, e.g., |T31 IT| [Til ITUj 
for the theory and applications of hyperbolic matrix and operator polynomials. 
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In this note we prove the following theorem. It states that the roots of combina- 
tions of hyperbolic matrix polynomials can be recast as eigenvalues of combinations of 
real symmetric matrices, under certain hypotheses. We denote by M the field of real 
numbers. 

Theorem 1.1 Let 

L{z) = Y^ LjZ^, M{z) = Y^ Mjz\ Me = Li = I, 

j=0 j=0 

be two monic n x n matrix polynomials such that 



aL{z) + (1 — a)M{z) is weakly hyperbolic for every a G 



(1.2) 



Assum,e in addition that the n x n matrix Li^i — Mf_i has n real eigenvalues (counted 
with multiplicities) . Then there exist ni x ni real symmetric matrices A and B such 
that for every a &M., the roots o/det (aL(z) + (1 — a)M{z)), counted according to their 
multiplicities, coincide with the eigenvalues of aA + (1 — cx)B, also counted according 
to their multiplicities. 

Conversely, if the roots o/det {aL{z) + (1 — a)M{z)), coincide with the eigenvalues 
of aA+{l —a)B {counted with multiplicities) for every a e M, where A and B are fixed 
real symmetric ni x ni matrices, then (ji.j^ holds and all eigenvalues of L£_i — Mf_i 
are real. 
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be the companion matrices of L{z) and of M{z), respectively. Then aCL + {l — Ci)CM is 
the companion matrix of aL{z) + (1 — a)M{z), and therefore the roots of det {aL{z) + 
(1 — a)M{z)) (counted with multiplicities) coincide with the eigenvalues of aCi + (1 — 
<^)Cm (also counted with multiplicities), for every a G M. (This is a standard fact in 
the theory of matrix polynomials, see, for example, j3].) Thus, aCL + (1 — C()Cm has 



n£ real eigenvalues for every a G M. Consider the homogeneous polynomial of three 
real variables a, /3, 7: 

P{a, /?, 7) := det (aC^ + ^Cm - ihi)- 

\i a + f3 7^ 0, the polynomial P(a, /3, 7) (as a polynomial of 7) has nl real roots (counted 
with multiplicities). If a + /3 = 0, then 

P(a, /3, 7) = ±7"^'-'^ ■ det (aL,_i + /3M,_i - 74) 

also has nl real roots, by hypothesis of the theorem. Thus, P(a,/3, 7) is hyperbolic in 
the direction of (0, 0, 1), in the sense of the Lax conjecture, see [HI HHj. By the main 
result of [13] (the proof in [THj is based on [HIEj), we have 

P(a, /3, 7) = det {aA + I3B - -fl) 

for some real symmetric matrices A and B. The direct statement of the theorem follows. 
To prove the converse statement, simply reverse the argument, taking into account 
that real symmetric matrices have all eigenvalues real. □ 

For further development of the theory of hyperbolic polynomials of several variables 
and many applications, in particular, mixed determinants, see [H]. 

2 Corollaries and applications 

We start by recalling Obreschkoff's theorem (see JHlEl)! which will be needed in the 
proof of the next corollary: Two real scalar polynomials f{z) and g{z) of degrees i and 
i — 1, respectively, have the property that f{z) + tg{z) has i real roots (counted with 
multiplicities) for every real t if and only if f{z) and g{z) have i and i — 1 real roots, 
respectively, and the roots of f{z) and of g{z) interlace (the cases when f{z) or g{z) 
have multiple roots and/or when f{z) and g{z) have common roots are not excluded 
here) . 

A proof of Obreschkoff's theorem can be given using the approach of Theorem 
II. H as follows (below, we formulate Obreschkoff's theorem in a slightly different but 
equivalent form): 

Proposition 2.1 The following statements are equivalent for scalar distinct monic 
relatively prime polynomials f{z) and h{z) of degree i: 

(1) The polynomials af{z) + l3h{z), a, [3 ^'^j a^ + /5^ 7^ , have all their roots real; 

(2) The polynomials otf{z) + (1 — a)h{z), a ^'^j have all their roots real; 



(3) f{z) has all its roots real, and the quotient h{z)/ f{z) has the form 

where Xj G M, and the real numbers Cj are all of the same sign; 

(4) Both f{z) and h{z) have i distinct real roots and the roots of f{z) and of h{z) 
interlace. 

Proof (1) clearly implies (2) . (3) and (4) are equivalent : Indeed, if (3) or (4) 
holds true, then f{z) has necessarily simple roots, and denoting the roots of f{z) by 
Ai < ■ ■ ■ < Ap, we see that in the representation ()2.ip . 



whereas 



..„...N_..„. M^) 


as 2; - 


-^ Aj, z > Xj, 


Gigni^Cjj sign 


.... .. ^ _ .... H^) 


as z - 


-> Aj+i, z < Aj+i 


sign {Cj^i) sign 



These equalities imply that that the roots of f{z) and h{z) interlace if and only if all 
the Cj's are of the same sign. 

We next prove that (2) implies (3). Arguing as in the proof of Theorem II. ![ we 
obtain that the characteristic polynomial oi aCf + f3Ch coincides with the characteristic 
polynomial of aA+f3B, for all a, /3 G M, where A and B are fixed (independent of a and 
(3) distinct real symmetric matrices. Taking a — (3 = we see that rank {A — B) < 1. 
Since polynomials / and h are distinct it follow that A = B ± xx^ for some nonzero 
vector X. Now 

f{z) = det {zl - Cf) = det {zl - A), 

h{z) = det {zl - Ch) = det {zl - 5), 

and 

jp- = det {{zl - B){zl - A)-^) = det (/ ± xx'^{zl - A)-^) = 1 ± x^{zl - A)-^x. 

This reduces, upon applying a diagonalizing real orthogonal transformation A 1— ;> 
U'^AU, and replacing x with U'^x, to ()2.H) with the real numbers Cj of the same 
sign, as required. 

Finally , let us prove prove the implication (3) => (1) . This means to prove 
that if (3) holds then for any real 7 the equation 44 + 7 = does not have roots with 
nonzero imaginary part . 



Concider a complex number z = a + bi , its real part Re{z) = a , its imaginery part 
Im{z) = b . li (3) holds then 



z - 



M^)^ME73T-'. 



where \j are real and the real numbers Cj are all of the same sign . Assume wlog that 
all Cj are positive . As 



M(.-A,)-) = ^^-^^-^^^. 



'h{z)^ 



thus we get that Im{j^) = -&Ej=i {a^x-)2+v^ ■ 

Therefore Im{j^) 7^ if Im{z) 7^ . This means that the equation j^ + 7 = does 

not have roots with nonzero imaginary part for all real 7. 

D 

We observe that checking condition (3) can be conveniently done using semidefinite 
programming. Indeed, let h{z) and f{z) be monic scalar polynomials with f{z) having 
all roots real, and consider a minimal realization 

^ = 1 + C{zl - Ar'B, 

fiz) 

where C, A, and B are real matrices. It is easy to see, using the uniqueness of a minimal 
realization up to a state isomorphism (similarity), that (3) holds, with the Cj's positive, 
if and only if there exists a positive definite matrix P such that 

AP = PA^, PC^ = B. 

The latter problem is a semidefinite programming problem. 

Another equivalent semidefinite programming problem is based on the following nice 

reformulation of Proposition 2.1 : 

The conditions of Proposition 2. 1 are equivalent to the existence of nonsingular real 

matrix D such that both DCfD~^ and DChD~^ are real symmetric . 

(Our proof of this statement is essentially the same rank one pertubration argument as 

in the proof of Proposition 2.1 ) . 

This gives the following semidefinite programming problem : 
is there exists a real positive definite P :^ such that 

PCf = CjP and PCh = ClP 

Our next corollary involves hyperbolic matrix polynomials. 



Corollary 2.2 Let L{z) he a hyperbolic matrix polynomial. Then there exist n£ x ni 
real symmetric matrices A and B such that the roots of det {L{z) +tL'{z)) coincide 
with the eigenvalues of A + tB {multiplicities counted), for every real number t. Here, 
L'{z) is the derivative of L{z) with respect to z. 

Proof. . By Obreschkoff's theorem, the matrix polynomial L{z) +tL'{z) is hyperbolic 
for every real t. Now apply Theorem 11.11 with M{z) = L{z) + L'{z). □ 

Note that the condition p.2j) implies (but is not equivalent to) the condition that 
every convex combination of L{z) and M{z) is weakly hyperbolic. It turns out that 
the latter condition can be conveniently expressed for hyperbolic matrix polynomials, 
which we will do next. 

Let L{X) be a hyperbolic n x n matrix polynomial. For every x G C", ||x|| = 1, let 

Ai(x)<A2(x)<--- ,<Xi{x) 
be the roots of equation ()1.1|) arranged in the nondecreasing order. The sets 

Aj{L) ■.= {Xj{x)\xeC'', \\x\\ =1}, 
called the spectral zones of -^(A), are obviously closed intervals on the real line: 

A,(L) = [57(L),5;(L)], j = l,2,...,£. 

A basic result in the theory of hyperbolic matrix and operator polynomials {^M Theo- 
rem 31.5], for example), states that two spectral zones either are disjoint, or have only 
one point in common. 

Proposition 2.3 Let L{X) and M{X) be two hyperbolic matrix polynomials of degree 
i. Then every convex combination aL{z) + (1 — a)M{z), < a < 1, is hyperbolic if 
and only if their spectral zones satisfy the inequalities 

max{5;(L),5;(M)} < mm{5T^,iL) , Sj^,{M)} , j = 1, ...,£- 1. 

For the proof apply Theorem 2.1]; this theorem gives necessary and sufficient 
conditions for all linear combinations of two given scalar polynomials to be hyperbolic. 

Using Theorem 11.11 and inequalities for eigenvalues of real symmetric matrices (see 
for example ^^), one can derive inequalities for eigenvalues of weakly hyperbolic matrix 
polynomials. We illustrate this for the case of the Horn inequalities. For a Hermitian 
m X m matrix X, we write its eigenvalues (repeated according to their multiplicities) 
in a non-decreasing order: 

Ai(X)<A2(X)<---<A^(X). 



An ordered triple {U, S, T) of nonempty subsets of {1,2,..., m} is said to be a Horn 
triple (with respect to m) if the cardinahties of [/, S, and T are the same, and the Horn 
inequalities 

hold true for every pair of Hermitian m x m matrices X and Y. A description of 
all Horn triples is known |HJ E] ; see also the surveys jlj |2] . For a weakly hyperbolic 
n X n matrix polynomial L{z) of degree i, we arrange the roots of det {L{z)) in the 
non- decreasing order: 

di{L)<d2iL)<...<dn,iL). 

Let T = {1 < ii < i2 < ■■■ < im ^ ni} C {1,2, ...,n£} , define T = {ni — i^ < 
ni — im-i < ■■■ < ni — ii} . 

Theorem 2.4 Let L{z) and M{z) he monic n x n matrix polynomials satisfying the 
hypotheses of Theorem \l.l\ Then for every Horn triple {U, S, T) with respect to ni, and 
for every a G M, the inequality 

J2 di {aL + (1 - a)M) < a i^ dj{L)\ + {I - a) { J] 4(M) j 

holds true. Here Sa = S if a >0 and S^ = S if a < 0. 

Proof. Let A and B be as in Theorem ll.il Then we have, using Theorem 11.11 and the 
Horn inequalities: 

^di{aL + {l-a)M)) = ^\i{aA+ {1 - a)B) 

ieu ieu 

< J]A,M) + ^A,((l-«)i?) 
jes keT 



a 



a 



J]A,(A)]+(l-a)[ J2 ^'^(B) 

Y^d,{L)\+{l-a)l Y, MM) 
i&Sc J \fceri_c 



and the proof is complete. 
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